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Abstract. Schubert calculus provides algebraic tools to solve enumerative prob 
lems. There have been several applied problems in systems theory, linear algebra 
and physics which were studied by means of Schubert calculus. The method is most 
powerful when the base field is algebraically closed. In this article we first review 
some of the successes Schubert calculus had in the past. Then we show how the 
problem of decoding of subspace codes used in random network coding can be for- 
mulated as a problem in Schubert calculus. Since for this application the base field 
has to be assumed to be a finite field new techniques will have to be developed in 
the future. 

1 Introduction 

Hermann Casar Hannibal Schubert (1848-1911) is considered the founder of enu 
merative geometry. He was a high school teacher in Hamburg, Germany. He studied 
questions of the type: Given four lines in projective three-space in general position, 
is there a line intersecting all given ones. This question can then be generalized to: 

Problem 1. Given N fc-dimensional subspaces % C C k+m . Is there a subspace 
"f c C k+m of complimentary dimension m = dimV such that 

f[]^i^{0},i=l,...,N. (1) 
Using a symbolic calculus he then came up with the following surprising result [21, 

79,]; 

Theorem 2. In case the subspaces % C C +m ,i = 1 , . . . ,N are in general position 
and in case N = km there exist exactly 

l!2!-(*-l)!(fen)! 
d( ' >~ m\(m+l)\---(m + k-iy: ( ) 

different m dimensional subspaces ~V C <C k+m which satisfy the intersection condi- 
tion (1). 

Note that two-dimensional subspaces in C 4 describe lines in projective space P 3 and 
Schubert hence claims in the case of four lines in three-space in general position that 
there are exactly d(2 1 2) = 2 lines intersecting all four given lines. 
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Schubert used in the derivation ot Theorem 2 Foncelet's principle ot preservation ot 
numbers which was not considered a theorem of mathematics at the time. For this 
reason Schubert's results were not accepted by the mathematics community of the 
19th century and Hilbert devoted the 15th of his famous 24 problems to the question 
if mathematicians can come up with rigorous techniques to prove or disprove the 
claims of Dr. Schubert. A rigorous verification of Theorem 2 was derived in the 
last century and we refer the interested reader to the survey article [14] by Kleiman, 
where the progress over time about Schubert calculus and the Hilbert problem 15 is 
described. 

In the sequel we introduce the most important concepts from Schubert calculus. 
Let F be an arbitrary field. Denote by Grass (k,n) = Grass (k,¥ n ) the Grassmann 
variety consisting of all ^-dimensional subspaces of the vector space F". Grass(£,n) 
can be embedded into projective space using the Plucker embedding: 

<p : Grass (k,F n ) — >.p(*H 
span(M[, . . . ,Uk) 1 — > F(«i A ... Am*). 

If one chooses a basis {e\,...,e n } of F" and the corresponding canonical basis of 
\ k F" 

{e h A...Ae ik \ I <h < ... <i k <n} 

then one has an induced map of the coordinates. If U is a k x n matrix whose row 
space rs(t/) describes the subspace % := span(Mi , . . . , u k ) and U[i\,..., i k \ denotes 
the submatrix of U given by the columns i\ , . . . , i k , then one readily verifies that the 
Plucker embedding is given in terms of coordinates via: 

rs(£7) i — > [det(f/[l,...,Jt]) : det(J7[l,...,Jfc- l,Jfc+l]) : ... : det(U[n -k+ 1, ...,«]). 

The k x k minors det(t/ [i \ , . . . , i*] ) of the matrix U are called the Plucker coordinates 
of the subspace % . 

The image of this embedding describes indeed a variety and the defining equations 
are given by the so called "shuffle relations" (see e.g. [15, 19]). The shuffle relations 
are a set of quadratic equations in terms of the Plucker coordinates. 
A flag & is a sequence of nested linear subspaces 

& : {0} C Vi C V 2 C . . . C V„ = ¥" 

tiaving the property that dim Vj = j for j = l,...,n. 
Denote by v = (vi , . . . , V^) an ordered index set satisfying 

1 < vi < . . . < Vfc < n. 

For every flag & one defines a Schubert variety 

5(v;^):={WeGrass(m,F") |dim(Wf|y v ,)>/ for i = 1, . . . ,k}. (3) 

The Schubert varieties are sub-varieties of the Grassmannian Grass(fc,F") and they 
contain a Zariski dense affine subset called Schubert cell and defined as: 

r(v:.y) :={W<=.Sf(v;.y) | dim(Wp|V v . i) = /-1: for/ I L\ (4) 
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In terms of Pliicker coordinates the defining equations of the Schubert variety S( v; 3£ 
are given by the quadratic shuffle relations describing the Grassmann variety to- 
gether with a set of linear equations (see [15]). 

A fundamental question in Schubert calculus is now the following: 

Problem 3. Given two Schubert varieties S(v; JP) and 5(v; J^"). Describe as explic- 
itly as possible the intersection variety 

S(v;jT)nS(v;#). 

Schubert's Theorem 2 can actually also be formulated as an intersection problem ol 
Schubert varieties. For this note that 

{V G Grass(fc,F* +m ) {0}} (5) 

describes a Schubert variety with regard to some flag and the theorem then states 
that in the intersection of N Schubert varieties of above type one finds d (k, m) m 
dimensional subspaces as solutions in general. 

In the case of an algebraically closed field one has rather precise information about 
this intersection variety. Topologically the intersection variety turns out to be a 
anion of Schubert varieties of lower dimension and the multiplicities are governed 
by the Littlewood-Richardson rule [9]. When the field is not algebraically closed 
much less is known. There has been work done over the real numbers by Frank 
Sottile [25, 26]. Over general fields very little is known and we will show in this 
article that the decoding of subspace codes can be viewed as a Schubert calculus 
problem over some finite field. The following example illustrates the concepts. 

Example 4. As a base field we take F = F2 = {0, 1} the binary field. Consider the 
Grassmannian Grass2(2,F 4 ) representing all lines in projective three-space P 3 . We 
would like to study Schubert's question in this situation: Given four lines in three 
space, is there always a line intersecting all four given ones. Clearly there are many 
situations where the answer is affirmative, e.g. when the lines already intersect 
in some point. In general this is however not the case as we now demonstrate 
Consider the following four lines in P 3 represented as row spaces of the following 
four matrices: 



10 
10 



10 
1 



110 1 
111 



1110 
10 11 



We claim that there exists no line in projective three-space P 3 , i.e. no two-dimensiom.l 
subspace in Grass2(2,F 4 ) intersecting all four given subspaces non-trivially. 

Grass2(2,F 4 ) is embedded in P 5 via the Pliicker embedding. Denote by 

u u :=detU[i,j],l <i<j< 4 
the Pliicker coordinates of some subspace £ Grass2 (2 , F 4 ) . The four lines impose 
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the linear constraints: 

M3.4 = 0, 

"1,2 = 0, 

Ml,2 + Ml 4 + U23 + U2,4 + M3 ; 4 = 0, 

Ml,2 + Ml,3 + «l,4 + «2,3 + W3,4 = 0. 

The points in P 5 representing the image of Grass2(2,F 4 ) are described by one 
quadratic equation (shuffle relation): 

Ml 2 M 3,4 + Ml 3M2 4 + Ml 4M2 3=0. 

Solving the 5 equations in the 6 unknowns results in one quadratic equation: 

(Ml, 4 ) 2 + M1 j4 M 2 ,3 + (M 2 ,3) 2 = 

which has no solutions over F 2 in P 5 . Note that there are exactly d(2,2) = 2 solu- 
tions over the algebraic closure as predicted by Schubert. 

Readers who want to know more on the subject of Schubert calculus will find mate- 
rial in the survey article [15]. 

The paper is structured as follows: In Section 2 we present results which were de- 
rived by Schubert calculus. In Section 3 we introduce the main topic of this paper, 
namely subspace codes used in random network coding. In Section 4 we show that 
list decoding of random network codes is a problem of Schubert calculus over some 
finite field. 

2 Results in Systems Theory and Linear Algebra De- 
rived by Means of Schubert Calculus 

In the past Schubert calculus has been a very powerful tool for several problem areas 
in the applied sciences. In this section we review two such problem areas and we 
show to what extend Schubert calculus led to strong existence results and better 
understanding. 

The pole placement problem 

One of the most prominent problems in mathematical systems theory has been the 
pole placement problem. In the static situation the problem can be described as 
follows: Consider a discrete linear system 

x(t + l)=Ax(t)+Bu(t), y{t)=Cx{t) (6) 

described by matrices A,B,C having size n X n, n X m and pxn respectively. Con- 
sider a monic polynomial 

q>(s) :=s" +a n _ l s"- 1 + --- + ai s + a Q e¥[s} 
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of degree n having coefficients in the base held IK. In its simplest version the pole 
placement problem asks for the existence of a feedback law u(t) = Ky(t) such that 
the resulting closed loop system 

x{t+\) = {A+BKC)x{t) (7) 

has characteristic polynomial (p(s). 

At first glance this problem looks like a problem from matrix theory whose solu- 
tion can be derived by means of linear algebra. Surprisingly, the problem is highly 
nonlinear and closely related to Schubert's Problem 1. This geometric connection 
was first realized in an interesting paper by Brockett and Byrnes [2] who showed 
that over the complex numbers arbitrary pole placement is generically possible as 
soon as n < mp and in case that the McMillan degree n is equal to mp then there 
are exactly d(m,p) static feedback laws resulting in the closed loop characteristic 
polynomial <p(s). The interested reader will find more details in a survey article by 
Byrnes [3]. 

The geometric insight one gained from the Grassmannian point of view was also 
helpful for deriving pole placement results over other base field. Over the reals 
the most striking result was obtained by A. Wang in [ ] where it was shown that 
arbitrary pole placement is possible with real compensators as soon as n < mp. Over 
a finite field some preliminary results were obtained by Gorla and the first authoi 
in | ]. 

U. Helmke in collaboration with X. Wang and the first author have been studying the 
pole placement problem in the situation when symmetries are involved [ I ()]. This 
problem then leads to a Schubert type problem in the Lagrangian Grassmannian. 

Sums of Hermitian matrices 

Given Hermitian matrices A\ , . . . ,A r G C" x " each with a fixed spectrum 

Ai(A/) > ... > A„(A/), 1 = 1,. ..,r (8) 

and arbitrary else. Is it possible to find then linear inequalities which describe the 
possible spectrum of the Hermitian matrix 

A r+ i :=Ai + ---+A r ? 

Questions of this type have a long history in operator theory and linear algebra. For 
example H. Weyl derived in 1912 the following famous inequality for any set of 
indices 1 < i, j < n with 1 < i + j— 1 < n: 

Aj+y-i(Ai +A 2 ) < Xi{M)+Xj{A 2 ). (9) 

In collaboration with U. Helmke the first author extended work by Johnson [ 1 3] and 
Thompson [28, 27] to derive a large set of eigenvalue inequalities. This was achieved 
through the use of Schubert calculus and we will say more in a moment. The ob- 
tained inequalities included in special cases not only the inequalities by H. Weyl but 
also the more extensive inequalities from I.idskii and Freede Thompson [28]. 
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In order to make the connection to Schubert calculus we follow [ J and denote 
with V\i,...,v n i the set of orthogonal eigenvectors of the Hermitian operator A/ 
! = l,...,r+l. 

Using these ordered set of eigenvectors one constructs for each Hermitian matrix A/ 
the flag: 

&l • {0}CV U CV2lC...CV n i = C n (10) 
defined through the property: 

Vml :=span(v 1/; ...,v m/ ) m = l,...,n. (11) 

The connection to Schubert calculus is now established by the following result as it 
can be found in [9]. The theorem generalizes earlier results by Freede and Thomp 
son [28]. 

Theorem 5. Let A\, . . . ,A,- be complex Hermitian n x n matrices and denote with 
&\ , . . . , & r + 1 the corresponding flags of eigenspaces defined by (11). Assume A r+ \ 
\\-\ YA r . and let l t = (in, . . . be r + 1 sequences of integers satisfying 

\<hi<...<iid<n, l = l,...,r+l. (12) 

Suppose the intersection of the r+ 1 Schubert subvarieties of Grass(&, C") is nonempy, 
'.e.: 

sfej^n.-n^wi^+i)^ - d3) 

Then the following matrix eigenvalue inequalities hold: 

k r k 

£ K-i Jr+l+ i(Ai + ■■■ +A r ) X 'A A i) (14) 

j=i i=i j=i 

k r k 

£ A,-,. +1 ( A i H +A r ) < t V(,+i (Ai). (15) 

j=l '=1.7=1 

In 1998 Klyachko could show that the inequalities coming from Schubert calculus 
as described in Theorem 5 are not only necessary but that they describe a Polytope 
of all possible inequalities. The interested reader will find Klyachko's result as well 
as much more in the survey article by Fulton [6]. 

A priori classical Schubert calculus provides very strong existence results. It is a 
different matter to derive effective numerical algorithms to compute the subspaces 
which satisfy the different Schubert conditions. For this reason Huber, Sottile and 
Sturmfels [12] developed effective numerical algorithms over the reals. As we will 
demonstrate in the next sections it would be very desirable to have effective numer 
ical algorithms also in the case of Schubert type problems defined over some finite 
field. 
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3 Random Network Coding 

In network coding one is looking at the transmission of information through a net- 
work with possibly several senders and several receivers. A lot of real-life applica- 
tions of network coding can be found, e.g. data streaming over the Internet, where a 
source wants to send the same information to many receivers at the same time. 

The network channel is represented by a directed graph with three different types of 
vertices, namely sources, i.e. vertices with no incoming edges, sinks, i.e. vertices 
with no outgoing edges, and inner nodes, i.e. vertices with incoming and outgo- 
ing edges. One assumes that at least one source and one sink exist. Under linear 
network coding the inner nodes are allowed to forward linear combinations of the 
incoming information vectors. The use of linear network coding possibly improves 
the transmission rate in comparison to just forwarding information at the inner nodes 
[1]. This can be illustrated in the example of the butterfly network: The source S 




Figure 1 : The butterfly network under the forwarding and the network coding model. 

wants to send the same information, a and b, to both receivers Rl and R2. Under 
forwarding every inner node forwards the incoming information and thus has to de- 
cide on either a or b (in this example on a) at the bottleneck vertex, marked above 
by x. Thus, Rl does not receive b. With linear network coding we allow the bottle- 
neck vertex to send the sum of the two incoming informations, which allows both 
receivers to recover both a and b with a simple operation. 

In this linear network coding setting, when the topology of the underlying network 
is unknown or time-varying, one speaks of random (linear) network coding. This 
setting was first studied in [11] and a mathematical model was introduced in [17], 
where the authors showed that it makes sense to use vector spaces instead of vectors 
over a finite field ¥ q as codewords. In this model the source injects a basis of the 
respective codeword into the network and the inner nodes forward a random linear 
combination of their incoming vectors. Therefore, each sink receives a linear combi- 
nations of the original vectors, which span the same vector space as the sent vectors, 
if no errors occurred during transmission. 

In coding practice the base field is a finite field F t/ having q elements, where q is 
a prime power. F* := F f/ \ {0} will denote the set of all invertible elements of V q . 
We will call the set of all subspaces of F^ the projective geometry of ¥ n q , denoted by 
&(q,n), and denote the Grassmannian Grass(A:,F^) by Grass^ (&:,«). 

There are two types of errors that may occur during transmission, a decrease in 
dimension which is called an erasure and an increase in dimension, called an inser- 
tion. Assume ^ G 3^{q,n) was sent and erasures and insertions occurred during 
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transmission, then the received word is of the type 

where ^ is a subspace of and S G ^(q,n) is the error space. A random net- 
work coding channel in which both insertions and erasures can happen is called an 
operator channel. 

In order to have a notion of decoding capability of some code a good metric is 
required on the set ,<3?(q,n): The subspace distance is a metric on ^(q^n) given by 

ds(W,y) = dim(^ + r)-dim(^nr) 

= dim(f/) +dim(r) - 2dim(^ n Y) 

for any ' ,"V G &(q,n). Another metric on £P(q,n) is the injection distance, de- 
fined as 

di{W,-y ) =max{dim(^),dim(r)} - dim(^ n f). 

Note, that for , f G Grass 9 (k, n) it holds that ds{W, V) = 2di(W, "¥) . A subspace 
code is simply a subset of £P(q,n). If ^ C GrasSq(k,n), we call it a constant 
dimension code. The minimum distance of a subspace code is defined in the usual 
way. 

Different constructions of subspace codes have been studied, e.g. in [4, 5, 16, 17, 
18, 20, 24, 30]. Some facts on isometry classes and automorphisms of these codes 
can be found in [29]. 

The set of all invertible n x n-matrices with entries in ¥ q , called the general linear 
group, is denoted by GL„. Moreover, the set of all k x n-matrices over F 9 is denoted 

by F* x ". 

Let U G F*j x " be a matrix of rank k and 

°l£ =rs(U) :=row space(C/) G Grass c/ (£,«). 

One can notice that the row space is invariant under GZ^-multiplication from the 
left, i.e. for any T G GL^ 

% = rs(U) =ts(TU). 

Thus, there are several matrices that represent a given subspace. A unique represen- 
tative of these matrices is the one in reduced row echelon form. Any k x n-matrix 
can be transformed into reduced row echelon form by a T G GL^. 

Given U G ¥ k q xn of rank k, G Grass 9 (/t,n) its row space and A G GL n , we define 

<%A :=rs(E/A). 

Let U,V G ¥ k q x " be matrices such that rs(U) = rs(V). Then one readily verifies than 
rs(t/A) = ts(VA) for any A G GL„, hence the operation is well defined. 
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Decoding subspace codes 

Given a subspace code ft C 3?>(q,n) and a received codeword ffl £ £P(q,n), a max- 
imum likelihood decoder decodes to a codeword % £ that maximizes the proba- 
bility 

P(M received | ^ sent) 

over all ^ £ ^. 

A minimum distance decoder chooses the closest codeword to the received word 
with respect to the subspace or injection distance. Let us assume that both the era- 
sure and the insertion probability is less than some fixed £. Then over an operatoi 
channel where the insertion probability is equal to the erasure probability, maximum 
likelihood decoding is equivalent to minimum distance decoding with respect to the 
subspace distance while in an adversarial model it is equivalent to minimum distance 
decoding with respect to the injection distance [23]. 

Assume the minimum (injection) distance of is d, then if there exists % £ with 
di(M, < ^y~, then °?/ is the unique closest codeword and the minimum distance 
decoder will always decode to °i/ . 

Note, that a minimum subspace distance decoder is equivalent to a minimum injec- 
tion distance decoder when c € is a constant dimension code. Since we will investi- 
gate constant dimension codes in the remainder of this paper we will always use the 
injection distance. All results can then be carried over to the subspace distance. 
A very important concept in coding theory is the problem of list decoding. It is the 
goal of list decoding to come up with an algorithm which allows one to compute all 
code words which are within some distance of some received subspace. 
For some % £ £?{q : n) we denote the ball of radius e with center in &(q,n) 
by B e (^/ ). If we want to describe the same ball inside Grass g (fc,n) we denote it by 
B k e (^ ). Note that for a constant dimension code the ball B k e (^ ) is nothing else than 
some Schubert variety of Gmss c/ (k,n). 

Given a subspace code c € C &(q,n) and a received codeword S% £ &(q,n), a list 
decoder with error bound e outputs a list of codewords ^ , . . . , % n £ c € whose in- 
jection (resp. subspace) distance from M is at most e. In other words, the list is 
equal to the set 

B e {M)r\^. 

If ^ is a constant dimension code, then the output of the list decoder becomes 

4 List Decoding in Plucker Coordinates 

As already mentioned before the balls of radius t (with respect to the subspace dis- 
tance) around some % £ Grass t/ (A:,n) forms a Schubert variety over a finite field. In 
terms of Plucker coordinates it is possible to give explicit equations. For it we need 
the following monomial order: 

(h,...,ik)>Uu-Jk) 3N£N> :ii=jiVl<Nandi N >j N . 
It is easy to compute the balls in the following special case. 
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Proposition 6. Define := rs[ 4 X /t ®kxn-k ]• Then 

B*(%) = {Y e Grass^n) \ <p{f) = [mi k ■ ••• : «], 

= V(ii,...,i' fc ) ^ (t + l,...,k,n-t+l,...,n)} 

Proof. For ^ to be inside the ball it has to hold that 

di(%,Y)<t 
k - dim{%) Di / )<t 
=> dim(^ nf)>k-t 

i.e. k — t many of the unit vectors e\ , have to be elements of V . Therefore 
(p(-f) hastofulfiU 

=0if (h,...,ik) i (t + l,...,k,n-t+l,...,n). 

a 

With the knowledge of B k t (ftt§) we can also express B k t (°?/ ) for any °^ G Grass 9 (£,«) 
For this note, that for any G Grass^ (k, n) there exists an A G GL n such that %A - 
% ' . Moreover, 

b!(%a)=b'[(%)a. 

For simplifying the computations we define (p on GL„, where we denote by A, t . ^ 
the submatrix of A that consists of the rows z'i, . . . , 4 

<p : GL„ — > ^^(1) 

detA^ ..^[1,...,^] ... detA[ k [n — k+l, 

detA n _ k+l ^ n [l,...,k] ... detA„_ jfe+ i ) .., ) „[n-&+l,...,rt] 
Lemma 7. Lef a i/ G Grass £/ (fc,n) a/iofA G GL n . It holds that 

<p(WA) = <p{<fr)(p(A). 
Theorem 8. Let = <%A G Grass^(/fc,«). Then 

B\{W)=B k t {%A) 

= {Y e Gmss q (k,n) | <p( K 1 / )<p(A^ 1 ) = [Mo,...,fc-i : • • • : Mn-*+i,...,»L 

Mi lr ..A = 0V(ii,...,4) £ (f+l,...,^,n-f+l,...,n)}. 

There are always several choices for A G GL„ such that ^oA = 9/ . Since GL^ is 
very large we try to choose A as simple as possible. We will now explain one such 
construction: 



1 . The first k rows of A are equal to the matrix representation II of a // 
10 
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2. Find the pivot columns or U (assume that U is in reduced row echelon form). 

3. Fill up the respective columns of A with zeros in the lower n — k rows. 

4. Fill up the remaining submatrix of size n — kxn — k with an identity matrix 

Then the inverse of A can be computed as follows: 

1 . Find a permutation (J £ S„ that permutes the columns of A such that 



(7(A) = 



2. Then the inverse of that matrix is 



(7(A)- 



/, U" 

I n -k 



h -u" 

I n -k 



3. Apply a on the rows of (7(A) The result is A l . One can easily see this if 
one represents a by amatrixS. Then one gets (SA)~ l S = A~ l S~ l S = A -1 

Example 9. In ^(2,4) we want to find 

B?(f) = {fe%(2,4) | fn^ = 1} 

for 

^ =rs(f/) =rs 
We find the pivot columns U [1 , 3] and build 



10 
11 



A = 



/ 1 \ 

11 

10 

\ 1 J 



Then we find the column permutation a = (23) such that 

(7(A) = 



/ 1 \ 

10 1 

10 

\ o o o i J 



Now we can easily invert as described above and see that (7(A) 1 
(7 on the rows and get 

^ 1 ^ 
10 
10 1 
\ 1 / 



(7(A). We apply 
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Then 































V 

































From Theorem 8 we know that 



B\{W) = S? 2 (2,4) I (ptfMA- 1 ) = [ M i,2 : ■ ■ ■ : Ma], 
Mi I ,i 2 =0V(ji,r 2 )^(2,4)} 
= {f G # 2 (2,4) I (p(y)(p(A- 1 ) = [ M i,2 : Mi,3 : - : M3, 4 ],M3,4 = 0} 



Now let (p{~f) = [vj.2 : V1.3 : Vi, 4 : v 2 , 3 : v 2 ,4 : v 3 , 4 ], then 

<p(r)<p(A _1 ) = [Vi )3 : Vi, 2 : Vi, 3 + v M : v 2 , 3 : v 3:4 : v 2 , 3 + v 2)4 ] 

and hence 



= {V G % (2,4) I p(r) = [vi, 2 : Vi, 3 : Vi, 4 : v 2 , 3 : v 2 a : v 3 , 4 ], V 2 , 3 + v 2 , 4 = 0} 
= { ^ G #2(2,4) I p(r) = [Vi, 2 : v L3 : Vi, 4 : v 2 , 3 : v 2 , 4 : v 3)4 ], v 2;3 = v 2;4 }- 

Note, that we do not have to compute the whole matrix <p(A~ l ) since in this case we 
only need the last column of it to find the equations that define B\(fy). 

5 Conclusion 

The article explains the importance of Schubert calculus in various areas of systems 
theory and linear algebra. The strongest results in Schubert calculus require that the 
base field is algebraically closed. The problem of list decoding subspace codes is a 
problem of Schubert calculus where the underlying field is a finite field. It will be a 
topic of future research to come up with efficient algorithms to tackle this problem 
computationally. 

Many of the results we describe in this paper were derived by the first author in col- 
laboration with Uwe Helmke. This collaboration was always very stimulating and 
the first author would like to thank Uwe Helmke for this continuing collaboration. 
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